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FOREWORD 


An  investigation  was  initiated  in  1950  of  the  problems  involved  in 
the  design  of  an  automatic  pressure  control  system  for  the  Engine  Test 
Facility  and  Ram- Jet  Addition,  Arnold  Engineering  Development  Center. 
An  analysis  of  the  plant  was  made,  and  the  performance  of  all  compo¬ 
nents  of  the  control  system  was  specified. 

Individual  responsibilities  for  different  phases  of  the  investigation 
are  as  follows: 

G.  V.  Schwent  -  determined  the  basic  schemes  of  analysis  and 
directed  the  study. 

D.  W.  Russell  -  supervised  work  on  distributed  parameter  equa¬ 
tions  and  completed  the  lumped  parameter  analysis. 

W.  K.  McGregor  and  L.  F.  Burns  -  derived  distributed  parameter 
equations  and  determined  the  frequency  response  solutions. 

Roger  Messick  -  obtained  the  indicial  solution  to  the  distributed 
parameter  equations. 

C.  T.  Coffee  -  contributed  suggestions  and  directed  the  mathema¬ 
tical  calculations. 

The  report  was  compiled  by  Mr.  McGregor. 
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ABSTRACT 


Methods  are  developed  for  the  prediction  of  frequency  and  indielal 
response  of  gas  flow  systems  on  both  a  lumped  parameter  and  a  dis¬ 
tributed  parameter  basis.  The  laws  of  physics  which  apply  to  gas  flow 
are  derived,  and  from  them  the  differential  equations  which  govern  the 
system  response  are  developed.  The  methods  of  operational  calculus 
are  used  to  solve  these  equations  for  the  case  of  sinusoidal  and  indicial 
forcings.  The  various  implications  of  these  solutions  with  their  bearing 
on  control  of  pressure  variables  are  discussed.  Examples  are  pre¬ 
sented  in  the  appendices  to  add  clarification  to  the  method  and  its  im¬ 
plications.  The  report  correlates  existing  knowledge  of  gas  flow 
systems  as  a  basis  for  formulating  an  inclusive  reference  to  be  used  in 
the  design  of  gas  flow  systems  having  desirable  dynamics  and  in  the 
design  of  control  systems. 
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NOMENCLATURE 


The  gravitational,  dimensional  system  will  be  used  throughout.  The  numeri¬ 
cal  unit  system  of  Appendix  3  will  use  the  British  units;  pounds  force,  feet, 
second,  degrees  Rankine. 


A 
A o 

A*,  B*,  C*,  D*.  E*,  F*  - 

C 

G 

K 

L 

M 

P 

Q 

R 

R* 

U 

V 

w 
w * 

7 


Area  of  valve  restriction 
Duct  cross-sectional  area 
Boundary  value  functions  of  s  or  jc u 
Lumped  pneumatic  capacitance 
Magnitude  ratio  (or  gain) 

Arbitrary  constant 
Total  length  of  ducting 
Mach  Number 
Absolute  pressure 
A  defined  function  of  s,  Q  =  (T  + 

\  hs 

Ideal  gas  constant;  In  Appeiidix  3,  R  -  53.3  fi/°R 

Residue 

Velocity  of  moving  fluid  in  a  duct 
Volume 

Weight  flow  of  gas 
Weight  of  gas 

Terminal  pneumatic  impedance 


a  Incremental  change  (or  small  perturbation)  in  valve 

area  from  the  steady  state  value,  A 

b  Used  in  definition  of  e 

c  Distributed  pneumatic  capacitance 
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f 

h 

i 

k 

n 

P 

r 

s 

t 

u 

U) 

X 

y 

z 

a,t 3 

y 

s 

e 

< 

7 

0 

e 

A 


Function  of  one  or  more  variables 
Distributed  pneumatic  inertance 
The  quantity  \j~-T 

Assigned  constant  denoted  by  subscript;  integer  without 
subscript 

Integer 

Incremental  change  (or  small  perturbation)  in  pressure 
from  the  steady  state  value,  P 

Distributed  pneumatic  resistance 

The  Laplace  complex  operator 

Time 

Velocity  of  pressure  wave  propagation  (velocity  of  sound) 

Incremental  change  (or  small  perturbation)  in  weight  flow 
from  the  steady  value,  W 

Distance 

Distributed  pneumatic  admittance 

Distributed  pneumatic  impedance 

The  real  numbers  defined  by  \fzy~  =  a  +  // 3 

Ratio  of  specific  heats  (y  =  1.  4  for  diatomic  gas) 

A  defined  quantity 

The  quantity  Lim  (1  +  x)X^x  ,  or  2.  178  ... 
x-*  o 

The  dimensionless  damping  ratio  of  a  second  order 
differential  equation 

Integer 

Absolute  temperature 

Incremental  change  (or  small  perturbation)  in  absolute 
temperature  about  the  steady  state  value  0 

Wave  length  of  a  sinusoidal  pressure  wave 
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n 

P 

a 

T 


6) 


The  factor  3.  1416  .  .  . 

Gas  density 

The  quantity  defined  by  e2an  =  — 
Time  constant 
Phase  angle 

Angular  frequency  of  oscillation 


z  -4f  o 

z  Q 


Subscripts: 

Numerical  subscripts  1,  2,  3  etc.,  and  alphabetical  subscripts  a,  b,  c,  etc., 
are  used  principally  to  distinguish  between  like  expressions  such  as  /,  k,  and  w. 

Alphabetical  subscripts  £,  o,  and  *  are  specifically  used  to  designate  stations 
--£  being  used  to  designate  inlet  and  o  to  designate  outlet. 

dp  * 

Dot  Notation  is  used  to  express  derivative  with  respect  to  time,  e.  g. ,  — —  =  P. 
-  at 
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INTRODUCTION 

Many  industrial  organizations  have  undertaken  extensive  analyses  of  their  pro¬ 
cesses  in  order  to  determine  more  exactly  their  control  requirements.  One  of 
these  processes  concerns  the  control  of  pressure  and  flow  of  gases  through  sys¬ 
tems  of  ducting,  referred  to  as  gas  flow  systems  in  this  paper.  An  example  of 
such  a  system  is  a  wind  tunnel  in  which  airframe  models  or  power  plants  are 
tested  under  accurately  simulated  Mach  number  and  altitude  conditions. 

Fundamental  to  the  design  of  a.  physical  plant  and  its  control  system  is  knowl¬ 
edge  of  the  behavior  of  variables  in  the  system  in  response  to  various  disturb¬ 
ances.  An  understanding  of  the  change  in  behavior  which  results  from  varying 
the  different  factors  and  properties  of  the  physical  design  is  also  required. 
Methods  of  acquiring  this  information  for  gas  flow  systems  are  developed  in  sub¬ 
sequent  sections. 

This  report  assembles  the  applicable  physical  relations  and  formulates  the 

response  of  pressure  in  gas  flow  systems  to  expected  disturbances.  There  does 

not  exist  a  treatment  of  this  problem  of  sufficient  scope  and  clarity  to  be  useful 

* 

to  the  control  designer.  To  be  sure,  the  laws  of  physics  which  apply  are  not  new 
and  may  be  found  in  various  textbooks  on  thermodynamics,  aerodynamics,  acous¬ 
tics  and  fluid  flow.  However,  the  objective  of  this  paper  is  the  use  of  these  laws 
and  their  formulation  into  a  useful  tool.  The  principal  applications  of  previous 
investigations  of  gas  flow  dynamics  have  been  concerned  with  such  things  as 
pressure  measurement  instrumentation  and  acoustic  muffler  design.  The  mathe¬ 
matical  similarity  between  the  laws  of  gas  flow  through  ducts  and  the  laws  of 
electricity  affords  an  analogy  that  has  been  used  advantageously  by  previous  in¬ 
vestigators  (Ref.  1). 
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The  principal  difference  between  this  work  and  existing  ones  on  the  same 
general  subject  is  that  here  the  gas  flow  rather  than  pressure  is  considered  as 
a  forcing  function.  This  is  necessary  in  an  analysis  for  control  system  design 
because  both  the  major  disturbance  and  the  manipulated  variable  are  gas  flows. 

The  methods  of  operational  calculus  are  used  in  this  report  because  of  the 
uniformity  of  expression  of  the  different  concepts  and  because  of  the  almost  uni¬ 
versal  use  of  the  Laplace  Transform  and  transfer  function  notation  in  the  synthe¬ 
sis  of  control  systems.  The  frequency  response  and  indicial  response  are  used  to 
express  the  behavior  of  a  system.  In  general,  the  treatment  here  utilizes  the 
methods  of  analysis  employed  by  texts  on  feedback  control  systems  and  thus  is 
directly  applicable  to  the  synthesis  problem. 
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THE  GOVERNING  LAWS  OF  PHYSICS 

Any  physical  system  is  governed  by  certain  laws  which  may  or  may  not  be 
formulated  in  mathematical  expression.  For  the  gas  flow  system  treated  herein 
these  laws  are  quite  simple  to  formulate  since  sufficient  assumptions  are  made 
so  that  "classical  mechanics"  hold  throughout.  By  a  gas  flow  system  is  meant  a 
connected  system  of  ducting  and  restrictions  through  which  a  gas  is  forced  to 
flow.  General  assumptions  which  must  be  made  are  as  follows: 

1.  The  flow  is  assumed  to  be  geometrically  one -dimensional  and  is  distrib¬ 
uted  continuously  and  uniformly  over  the  whole  duct  area,  being  a  con¬ 
tinuum. 

2.  Dissipation  by  radiation  and  thermal  conductivity  is  neglected. 

3.  The  fluid  is  a  pure  diatomic  gas  which  obeys  the  ideal  equation  of  state. 

The  pressure  of  a  gas  flowing  in  a  duct  is  a  function  of  both  time  and  dis¬ 
tance.  If  the  space  variable  is  neglected,  the  energy  storage  is  said  to  be 
"lumped",  and  the  system  is  said  to  be  a  "lumped  parameter"  system.  If  the 
space  variable  is  not  neglected,  the  system  is  a  "distributed  parameter"  system. 
The  equations  resulting  from  an  analysis  of  a  "lumped  parameter"  system  are 
ordinary  differential  equations  with  respect  to  time;  those  resulting  from  a 
"distributed  parameter"  analysis  are  partial  differential  equations  with  respect 
to  both  time  and  distance.  In  general,  the  coefficients  of  these  differential  equa¬ 
tions  are  not  constant  but  depend  on  the  pressure,  flow,  and  temperature 
in  the  system.  The  small  perturbation  theory  is  utilized  in  order  that  these 
coefficients  may  be  treated  as  constants  for  a  given  operating  point. 

THE  LUMPED  PARAMETER  EQUATIONS 

Three  elements  are  encountered  in  a  lumped  parameter  system:  energy 
storage  in  a  volumetric  capacitance,  flow  through  a  compressor,  and  flow  through 
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a  restriction.  The  relations  which  govern  flow  and  pressure  in  these  elements 
are  derived  as  follows: 

Pressure  -  Weight  Flow  Relationship  in  a  Volume 

Consider  the  volume  V  in  Fig.  1  through 
which  agas  at  pressure  P  and  temperature  0 
is  flowing  at  a  rate  K. 

The  assumption  is  made  that  tempera¬ 
ture  and  pressure  are  instantaneously  the 
same  anywhere  in  the  volume.  Then,  the 


W: 


P  0/3 


Volume  V 


W, 


FIG.  1.  SCHEMATIC  OF  VOLUME  V 


rate  of  change  of  density  of  the  gas  within  the  volume  is  the  difference  in  weight 
flow  into  and  out  of  the  volume  divided  by  the  volume.  That  is. 


do  Wi  ~ 

dt  V 


(1) 


The  equation  of  state  for  a  perfect  gas  is 


P  = 


RQ 


(2) 


Differentiating  equation  (2)  with  respect  to  time  and  combining  with  equation  (1) 
yields  the  general  relation  of  pressure  and  flow  in  a  volume. 


W,-  Vn  = 


RQ 


(P-gW  , 


(3) 


where  the  dot  notation  is  used  to  denote  derivative  with  respect  to  time.  If  the 
changes  of  state  of  the  gas  may  be  considered  as  isentropic  processes,  then 

y  (4) 


and  equation  (3)  reduces  to 


y  -  I 

p@  =  const  , 


W.  _  w  - _ _ _ 

1  yRQ 


P  , 


(5) 


This  latter  condition  that  the  processes  be  isentropic  can  be  satisfied  by  al¬ 
lowing  the  variables  to  change  in  small  perturbations  about  some  steady  state 
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value.  That  is, 

P  =  Pss  +P 
1 *i  =  Wss  +  wi 

W0  =  %ss  +  U)0 

®  =  ®ss  +  0  • 
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If  these  perturbations,  denoted  by  the  lower  case  letters,  are  sufficiently  small 
the  processes  are  essentially  isentropic.  Making  this  assumption  and  defining 


we  obtain 


C  = 


V 

yP®ss 


p  -L.(wi  -  w0).  (6) 

On  making  the  Laplace  transformation  on  the  variables  of  equation  (6), 

p(s)  =  —  f  wi  (s)  -  (s)| .  ^ 

Cs  L  J 

The  quantity  C  is  the  lumped,  pneumatic  capacitance  of  volume  V.  If  0  is  small 
this  capacitance  will  be  constant.  It  should  also  be  noted  that  if  the  process  were 
isothermal  rather  than  isentropic  the  ratio  of  specific  heats,  y,  would  be  re¬ 
placed  by  unity.  Thus,  C  may  vary  by  a  factor  from  unity  to  1.4  depending  on  the 
degree  to  which  the  process  can  be  said  to  be  isentropic. 


Pressure  -  Weight  Flow  Relationship  in  a  Compressor 


In  any  gas  flow  system  there  must  be  some  method  of  forcing  the  gas  to  flow 
through  the  volume.  The  device  which  accomplishes  this  is  usually  referred  to 
as  a  compressor  if  it  forces  gas  into  the  volume  and 
as  an  exhauster  if  it  draws  gas  out  of  the  volume. 

The  flow  through  this  device  is  a  function  of  many  ri,vi 
variables.  Each  compressor,  or  exhauster,  is  char¬ 
acterized  by  a  unique  set  of  performance  curves  which 

FIG.  2.  SCHEMATIC  OF  COMPRESSOR 

show  flow  as  a  function  of  the  pressure  ratio  across 
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it  for  different  values  of  inlet  temperature.  Thus, 

V-flPi,  Po^O- 

Taking  the  differential  of  flow  yields 


dPi  dP0  d&i 


(8) 


(9) 


Making  use  of  the  theory  of  small  perturbations  so  that  linearity  is  assured 


w  =  keP i  +  kcP2  +  kt  6  ,  (10) 

where  the  following  definitions  have  been  made: 


w  -  dW 

ke  - 

dW 

Wi 

Pi  =  dPt 

dW 

h 

*Po 

Po  =  dPo 

e  -  d®i 

Ic.  - 

dW 

Kt  “ 

s®{  * 

The  constants  denoted  by  k  must  then  be  evaluated  at  each  operating  condition 
from  performance  curves. 


Pressure  -  Weight  Flow  Relationship  in  a  Restriction 

Control  of  pressure  in  a  gas  flow  system  is  accomplished  by  control  of  flow 
at  some  point;  this  in  turn  is  accomplished  by  controlling  the  area  of  some  re¬ 
striction,  usually  a  valve  of  some  configuration,  at  that  point.  In  general,  the 
flow  through  a  restriction  is  dependent  on  - 

Pi*0i  -  . 

upstream  pressure,  downstream  pressure,  - '  W 

temperature,  and  the  open  area  of  the  restric-  fig.  3.  schematic  of  restriction 

tion.  That  is. 


w  =f(Pi,P0,®i,A) 


(ID 


Taking  the  differential  and  using  the  theory  of  small  perturbations  as  before, 

w  =  pi  +  k0p0  +  kg  6  +  &0<*,  (12) 

where  the  following  definitions  have  been  made: 
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tv  =  dW 

ft*  = 

dW 

dpi 

Pi  =  dPi 

K  = 

ill V 

dp0 

Po  =  dPo 

h  = 

dW 

d®i 

6  =  d®i 

ka  = 

dW 

dA 

a  =  dA  . 

These  constants  depend  upon  knowledge  of  the  valve's  pres  sure- flow  - 
temperature  relationship,  which  may  be  expressed  either  by  families  of  curves 
or  by  approximate  mathematical  expressions,  A  method  is  given  in  Appendix  1 
whereby  these  constants  may  be  determined  for  some  types  of  valves.  This  meth¬ 
od,  with  sufficient  experimental  data,  can  be  extended  to  any  valve  configuration. 


THE  DISTRIBUTED  PARAMETER  EQUATIONS 

Consider  an  element  dx  of  a  straight  duct  of  constant  cross-sectional  area  A0 
in  which  a  compressible  fluid  of  density  p  and  at  pressure  P  is  flowing  at  a  rate 
W  and  velocity  V  as  shown  in  Fig.  4.  The  quantities  are  shown  at  any  instant  of 
time  t.  The  partial  differential 


equations  which  describe  the  pres¬ 
sure-flow  relationship  are  derivable 
from  an  expression  of  the  continuity 
equation  and  a  net  force  equation. 


P 

P 


dx 


W 

u 


I 


P  +  ~  dx 
dx 

1 1 +  f? dx 

Tw  +  f^dx 


FIG.  4.  SCHEMATIC  OF  INCREMENT  OF  FLOW  IN  A  DUCT 
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The  continuity  equation  states  that  the  rate  of  decrease  of  weight  within  the 
element  of  volume  A0dx  is  equal  to  the  net  rate  of  weight  flow  in  the  element, 

♦£■*•>-»  •  (13> 
The  weight  of  gas  in  the  volume  is  the  density  times  the  volume  and  hence, 

W*  =  (p  +  ~§~~  dx)AQdx  , 

Neglecting  the  higher  order  terms  yields 


J¥1  m. JP-AJx 

at  at  0 


Equating  (14)  and  (13)  yields 


aw 


-  Aq 


dp 


(14) 


(15) 


dx  v  at 

The  process  is  assumed  isentropic  as  in  the  case  of  lumped  parameter 
equations.  For  isentropic  processes  all  fluid  state's  have  the  same  entropy. 


This  fact  can  be  stated  as  follows  (see  Ref.  3): 
dP  dP  U 1 


dp 


e 


(16) 


dp  /constant) 
l  entropy  J 

where  u  is  the  velocity  of  propagation  of  a  pressure  disturbance.  The  equation 


of  continuity  can  then  be  stated  as 
dW  A0g  sP 


dP 


dx  u2  Bt  Bt 


(17) 


The  net  force  acting  on  the  fluid  within  the  elemental  volume  is  equal  to  the 

mass  of  the  fluid  times  acceleration.  The  mass  of  fluid  in  the  volume  is 

—  =  J-(p  +  ^2—  dx)A0dx- 
g  g  dx 

The  velocity  of  the  incremental  volume  is 


V  = 


V 


10P 


a  ,  w  . , 

+  - -  (- - )dx 


Bx 


*oP 


and  the  corresponding  acceleration  is 


BU 
d  t 

The  net  force  is  then 


a  v  _a 

at  A0p  a 


H 


w 


dx 


AoP 


■)  dx 


]■ 


PAC  _[>  .  *  A0  J-i-lp  1  (-i-(-I-)  toll  (18) 

L  dx  J  Lg  dx  J  U{  A0P  dt  L  Bx  A0p  J  J 
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Upon  simplifying  equation  (18)  and  neglecting  the  higher  order  terms  in  dx , 
equation  (18)  reduces  to 

_dP_=_P_  4p)  _  1  _  JL_  (19) 

dx  Aog  St  gA0  St  A0pg  St 

In  order  to  handle  equation  (19)  advantageously,  it  is  necessary  to  assume  that 

total  velocity  of  the  gas  stream  is  small.  Then,  the  last  term  of  equation  (19) 

can  be  neglected  and 

_  SP  =  _1 _  SW_  >  (20) 

dx  A  Qg  St 

If  the  quantities  are  treated  as  small  perturbations  about  a  steady  state 
value  as  in  the  lumped  parameter  treatment,  then  equations  (17)  and  (20)  become 


S  w 

-  =  c 

dp 

(21) 

dx 

St 

dP  -  h 

Sw 

(22) 

dx 

St 

There  also  occurs  a  loss  in  pressure  with  distance  due  to  various  resis¬ 
tances.  To  include  this  term  apparently  contradicts  the  assumption  of  isentropic 
processes  made  previously.  However,  change  in  enthropy  with  distance  does  not 
preclude  an  isenthropic  compression  in  an  incremental  element.  This  pressure 
drop  along  the  duct  can  be  taken  as  being  a  function  of  flow.  Hence,  over  a.  small 
incremental  change  (again  using  the  small  perturbations): 

dP  _  dP  dW 

dx  dW  dx 

or 

-  -i£-.r»;  ■  <23> 

dx  \W!\x 

Combiningthis  additional  pressure  drop  with  that  of  equation  (22)  yields 

_  =  hJOL-  +  rw  .  (24) 

Sx  dt 

The  quantities  c,  h  and  r  are  usually  referred  to  as  the  distributed  pneumatic 
capacitance,  inertance  and  resistance,  respectively. 
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A  SIMPLE  SYSTEM 


The  use  of  the  foregoing  physical  laws  can  best  be  illustrated  by  means  of 
an  example.  It  is  in  the  response  of  a  system  that  the  correlation  between 
the  lumped  and  distributed  parameter  descriptions  of  system  dynamics  can  best 
be  seen.  The  treatment  will  also  be  useful  since  the  simple  system  to  be  used 
is  a  sufficient  description  of  many  processes. 

The  system  to  be  studied  is  shown  in  Fig.  5.  In  Fig.  5(a)  the  system  is 
represented  as  a  lumped  parameter  system.  A  gas  flows  into  volume  V  through 
some  restriction  or  device 


at  a  rate  IfA  and  out  through 
another  device  at  a  rate  If  , 
with  a  pressure  P  and  a 
temperature®.  In  Fig.  5(b) 
the  system  is  represented 
as  a  distributed  parameter 


p,e 


W; 


V 


W, 


FIG.  5<o).  LUMPED  PARAMETER  SYSTEM 


.  system.  The  gas  flows  into 
a  straight  duct,  of  Length  L 
and  of  cross-sectional  area 
40(such  that  V  =  A 0L),  from 


dx 


i 


W; 


X 


I  I 

FIG.  5(b).  DISTRIBUTED  PARAMETER  SYSTEM 


w , 


the  same  source  as  in  5(a)  and  out  through  the  same  device  as  in  5(a),  maintaining 
the  same  steady  state  pressure  and  temperature. 

The  flow  into  the  system  is  an  independent  variable.  The  flow  out  of  the 
system  is  considered  to  be  proportional  only  to  the  pressure  P.  Temperature 
is  considered  constant  in  the  evaluation  of  the  capacitance  term. 
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PRESSURE  RESPONSE  TO  AN  INFLOW  DISTURBANCE  AS  DETERMINED  BY  THE 
LUMPED  PARAMETER  EQUATIONS 


Two  equations  describe  this  system  according  to  the  linearized,  small  per¬ 
turbation  theory  presented  previously: 


p(s)  = 


1 

■Cs 


(s)  -  i«0(s) 


and 


(25) 


%(s)  =  k0p(s) » 


(26) 


where  C  =  -  and  k0  = 

yRQ 

system  transfer  function. 


W 

— — -  •  Combining  these  equations  yields  the 
“ss 


P-(s) 

w.(s) 


lAo 


c  ^ 

T~s  +  1 

Ko 


,rs  +  1 


The  gain  of  the  system  is  K 

Frequency  Response 


V 


Ji-  and  the  time  constant  is 


ss 


(27) 


K  yR®Wss 


Of  interest  in  the  synthesis  of  control  systems  is  the  frequency  response 
of  the  system  to  be  controlled.  Replacing  s  by  ja  the  frequency  response  of 
transfer  function  (27)  can  be  determined.  The  straight  line  approximation  of 
this  frequency  response  is  shown  in  non-dimensionalized  form  in  Fig.  6. 


FIG.  6.  LUMPED  PARAMETER  FREQUENCY  RESPONSE 
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Indicia!  Response 


The  usual  manner  of  showing  the  time  response  of  a  system  is  to  express 
the  time  solution  to  a  step  input,  sometimes  referred  to  as  the  indicial  response. 
In  mathematical  form  this  is 


(28) 


The  graph  is  shown  in  Fig.  7. 


FIG.  7.  LUMPED  PARAMETER  INDICIAL  RESPONSE 


PRESSURE  RESPONSE  TO  AN  INFLOW  DISTURBANCE  AS  DETERMINED  BY  THE 
DISTRIBUTED  PARAMETER  EQUATIONS 

The  solution  of  equations  (21)  and  (24)  to  sinusoidal  and  indicial  forcings 
yields  the  distributed  parameter  response  of  the  example  system.  These  two 
solutions  are  sufficiently  dissimilar  to  require  separate  treatments.  The  first 
uses  the  standard  jo  operator,  whereas  the  latter  must  employ  the  Laplace 
operator. 


Frequency  Response 

Suppose  the  quantities  w  and  p  of  equations  (21)  and  (24)  to  be  small  sinusoi¬ 


dal  perturbations.  The  impedance  notation  commonly  used  in  the  study  of  steady 

state  electrical  alternating  current  can  then  be  used  and  the  equations  become 

-,iw—  =  C  -  =  jcocp  =  yp  <29> 

dx  dt 


dp  du> 

-  =  tw  +  h  -= —  = 

dx  dt 

v/here  y  and  z  are  the  pneumatic  admittance 


(r  +  joh)w  =  zw  t  (30) 

and  impedance,  respectively,  while 
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joi  is  the  complex  frequency  operator.  It  should  be  noted  also  that  equations  (29) 
and  (30)  differ  from  equations  (21)  and  (24)  by  a  minus  sign  since  here  it  is  con¬ 
venient  to  measure  distance  from  the  outlet  end  of  the  system. 

The  analogy  between  the  distributed  parameter  gas  flow  system  and  an  elec¬ 
trical  transmission  line  can  be  observed  through  examination  of  equations  (29) 
and  (30). 

Taking  the  derivative  of  equation  (30)  with  respect  to  distance  yields 


d 2  p 
dx 2 


dw 

dx 


=  zyp 


or 


d2  p 
dx2 


-  zyp 


=  o 


A  solution  to  this  equation  is 


Px  =  A* e 


x 


+  B*  e 


~x\fzy 


(31) 


(32) 


(33) 


and  from  equation  (29) 

where  A*  and  B*  are  to  be  determined  from  boundary  conditions.  For  the  system 

under  consideration,  at  *  =  0  ,  px  =  p0  and  wx  =  w0 .  After  evaluating  A*  and  B* 

and  simplifying,  the  solution  can  be  written  as 

px  (/<*>)  =  pQ  cosh  x\J~zy  +  u>0  ■yjy-~  sinh  xyj  zy  (34) 

u>x  (jc»)  =  u>0  cosh  Xyfzy~  +  poyJ—  sinh  Xyfzy  . 


As  mentioned  previously  the  disturbance  to  the  system  is  a  flow  oscillation. 


U) 


i » 


at  *  =  L.  Hence,  we  desire  the  pressure  response  to  this  disturbance  at  any  po¬ 
sition  in  the  ducting.  Making  the  definitions 


Z  = 


JL 

dW„ 


O  _ 


Po 


’  o  0 

\[W  m  a  +  iP  » 


=  (0  -  ja)  . 

coc  ^  J  » 
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where 


V 


%  coc  [w  h  +V&Ji  h1  +  r1] 


arc 

a  =  - - 

2ft  * 

and  using  certain  trigonometric  and  hyperbolic  identities,  equations  (34)  and  (35) 


may  be  manipulated  to  give 


where 


Px  <&>  m  ft  +  j  f2 
wi  (jta)  ft  +  j  ft 


Gx 


+  u1 
+  A* 


(36) 

(37) 


4%  =  Tan"*  — h 
fi 

and  the  functions  of  ja  are  given  by 


Tan” 


ft 


(38) 


ft  =  Z  cos  (Bx  )  cosh  (ox)  +  — — —  cos  (px)  sinh  (ax)  +  — - —  sin  (Bx)  cosh  (ax) 

coc  coc 

f2  =  Z  sin  (Bx)  sinh  (ax)  +  — — —  sin  (Bx)  cosh  (ax)  -  — - —  cos  (Bx)  sinh  (ax) 

coc  coc 

/j  =  cos  (ftL)  cosh  (aL)  +  j  ft  ft2  cos  sinh  ( ° £) - a*  '+  p*  sin  (/9L)  cosh  (aL) 

ft  =  sin  (/3L)  sinh  (aL)  +  2  sin  (/3 L)  cosh  (aL)  + — f  cos  (|8L)  sinh  (aL)  . 

a  +  p  a  +  o 

Thus,  equations  (36),  (37)  and  (38)  give  the  frequency  response  of  the  system 
of  Fig.  5  (b)  at  any  position  x,  where  x  is  measured  from  the  outlet. 

Considerable  simplification  of  this  complicated  frequency  response  function 
results  for  the  special  case  of  zero  resistance  (r  =  0) .  Then  ft  =  co/u,  a  =  0  and 
the  functions  of  ja  become 

fx  =  Z  cos  fix  =  Z  cos  -St-  x  ft  ~  cos  fiC  =  cos  G 

f,  =  sin  ftx  =  —5—  sin  x  f.  =  Z  afj . sin  ftL  =  Z  cu  sin  -2L-  L 

12  ac  ^  uc  u  14  ft  u 

At  x  =  o,  ft  =  Z  and  f2  =  o.  The  functions  ft  and  ft  range  from  maxima  to 
minima  at  ftL  =  n  v/2  ,  where  n  is  any  integer  and  for  r  =  0,  ft  =  owing  to  the 
nature  of  the  trigonometric  functions.  This  occurs  at  such  values  of  frequency, 
a>,  such  that  L  is  equal  to  some  multiple  of  a  quarter  wave  length.  That  is,  at 
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L  =  -T-  X  -  -a— ,  and(Dc=  V”  -  where  <mc  is  denoted  as  a  "critical"  frequency. 
From  this  alone,  the  maxima  and  minima  of  the  frequency  response  can  be  deter¬ 
mined.  This  treatment  yields  the  same  information  when  applied  to  the  position 
x  =  L,  Thus,  the  frequency  response  plots  of  the  distributed  parameter  system 
can  be  approximately  determined  as  shown  in  Fig.  8. 

Indirial  Response 

The  methods  of  the  transformation  calculus  can  be  used  to  obtain  the  pressure 
response  to  a  step  disturbance  in  flow.  By  obtaining  the  Laplacian  transform  of 
response,  finding  the  singular  points  and  residues  of  the  transformed  equation 
and  applying  the  "Inversion  Theorem"  (see  Ref.  4),  the  time  solution  can  be  e- 
valuated.  This  process  will  be  applied  in  this  section  to  equations  (21)  and  (24), 
since  here  it  is  convenient  to  measure  *  from  the  inlet. 


Differentiating  equations  (21)  and  (24)  with  respect  to  distance  and  time  re¬ 
spectively  yields 


d2  w 

d2  P 

(39) 

dx2 

dt  dx  ‘ 

_  d>P  - 

dxdt 

,  d2w  dw 

°h  dt2  +r  dt  ; 

(40) 

or,  upon  combining 

d2w 
dx 1 

,  d2w  dui 

=  hc  dt2  +  r°  dt  • 

(41) 

In  a  similar  manner,  differentiating  equations  (21)  and  (24)  with  respect  to  time  and 

distance,  respectively,  and  combining, 

yields 

d2P 
dx 2 

"  hc  dt2  +  rc  dt  ’ 

(42) 

If,  initially. 

w(o,x) 

- pM  =  (ir)t  =  o  -  0i 

then  the  Laplace  transformation  of  equations  (41)  and  (42)  yields 

d2w  (s,x) 

=  (he  s  +  rc)  s  w(s,x)  , 

(43) 

dx 2 

d2  p(s,x) 

=  (hc  s  +  rc)  s  p(s,x )  . 

(44) 
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These  equations  have  the  solutions 

w(s,x)  =  C  *e  < +  ™>V’  *  +  D*e-(Acfi2  +  rcs),/l  * 
p(S,*)  =  E*e  (Acs‘  +  TCs)V’  x  +  F*e~{kcs2  +  X  . 


(45) 

(46) 


The  total  length  of  the  system  is  taken  as  L.  At  *  =  L,  w(s)  =  u>o  (s),  and 
p(s)  =  p0(s) .  In  addition,  the  "terminal  impedance"  is  defined  by  Z  =  , 

and  it  is  also  expedient  to  make  the  definition 

Q  ,  (jUjJ±\ 

\  hs  / 

Upon  evaluating  C*  ,  D* ,  E *  and  F*  as  functions  of  s,  equations  (45)  and  (46)  be¬ 


come 


w(s;x)  = 


u>0  (s) 


2  Qyff 


Q  Vhc  s(L  -  x) 


-  (z-Qyff)' 


-  Qy/hc  s(L-  x) 


(47) 


p(s,*)  = 


w0  (s) 


Q  V  he  s(L-  x)  + 


(z-Qyff), 


-Q  \JTc  s{L-  x) 


.  (48) 


The  transfer  function  of  pressure  to  flow  disturbance  can  then  be  written  as 


(z-e# 


-  Q  V  he  s  {L  -  x) 


^  Q 

(Z  +  Q^l)eQ^s^-xK  (z-Qyff)  e-Q^s(L-x) 

If  the  disturbance  occurs  at  x  =  o,  the  transfer  function  can  be  written  as 

fz  +  Q  \  Q  'fhc  s(2L-x)  Q  \fhc  s  x 


LM.  m  yj±  Q  YlA 

'• >{s,0 )  v  -  v 


+  e 


+  Q\[f\  2<?  VAc  SL 


It  is  convenient,  here,  to  introduce  the  relation 


2  an 


Z-  #  g 

Z*  %/f  « 


(49) 


(50) 
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Then 


p  (s^x)  _  -IT  q 

ni  c  A  v  r  ' 


Q\!hc  s(2Z/-;t)  -2  no 


Q  TT 


U)(s,  o) 


2 Q  sfhTsL  -  2  077 
e  “1 


(51) 


It  is  now  expedient,  if  the  mathematics  is  to  be  kept  simple,  to  examine  the 

case  r  =  o,  therein  making  Q  equal  to  unity.  The  solution  is  desired  to  a  step  in 

flow;  that  is,  ug  (s)  =  •  The  poles  of  equation  (51)  are  (1)  s  =  0  and 

s 

/o\  au  x  <.  Zs'fhc  L  -  2on  2tr  ft  ,  n 

(2)  the  roots  of  e  =  e  =1,  which  are  s  =  — — iL. -  ( o  +  jk ), 

V  he  L 

where  k  =  0,  ±1,  ±2,  — 


The  residues  at  these  singularities  are  then 
1.  R*  =  Z 


2-r'T4F 


■ag-  ^  +  jk) 


+  e 


-77  —  (o  +  jk) 

h 


277 

jknx 

=  —  a~)'k  e  L 
n  *  c  o'+k* 


a  +  jk 


cosh 


nax 


where  the  simplification  of  R*  is  obtained  by  using  2njkx/L  m  (e2njk)L  =  y 

The  inverse  of  the  transform  is  obtained  by  summing  the  product  R*est  at 
each  singularity.  Then 

k  -  +■  < 

^JR*eSt  „  Z  + 

k  TT 


p(x,  t) 


S=  +  oo 

V  ±JIjl=JL 

Zj  ”  c  °i  +  kl 


— — - {axjk)  t 

inkx/Lc08h  JL£l  e  .  (52) 


By  introducing  the  identify 


jnkx/L  tfk  -j-qnk  j  nkx/L 

e  =  (-1)  e  e 


ink 

k  L 
=  (-1)*  e 


(x-7)  L) 


equation  (52)  becomes 


not 


p(x,  t) 


=  Z  + 


n  '  c 


cosh 


nax 


\/hc  L 


k  =  +  ' 


g  ~  jk  ,  ■ 
°'+k' 


i, 1  =  1,  3,  5  -  -  -  . 

ikn  [{  +  \fhc~  (*  - 17  Z,)] 


l)fe  e  The  L 


(53) 


k  = 
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Expanding  equation  (53)  yields 


P(*»  f>  =  Z  +  ±  JI  cosh  e^L  J  2a  _J_  +  1)  * 

tt  ^  c  L  I  2a2  /  /  a2  + 


.  A;  &  sin  k  6  * 


where  0*  =  (*  -  r)L)  ]  .  This  form  brings  to  attention  two  Fourier 

V  he  L 

series  identities  which  can  be  written  as 
C08h  O0*  .  JF-sinh  OT[-^7  +  ^(-l)k  -?“*?"]  * 


sinh  a  9*  =  - 


■f  .i»h  „  [  2  w>‘  ■ 


where ~n<9*<n.  Substitution  of  these  identities  (Ref.  5)  into  equation  (54)  yields, 
after  simplification, 


p(*»  t) 


Z  +  \  JL  e 


- —  rr — w  cosh  naX 

i  e'/Aci  L  e-o9* 


The  boundary  placed  on  9*  requires  that  for  each  rj,  t  be  confined  to  the  region 


V  he  0?-l)  -  *  J  <t  <  Vhc"  (j/  +  1)  -*j 


Allowing  rf  =  2  n+1  where  n  =  0,  1,  2,  —  and  remembering  the  definition  of  2 an  for 


<2=1,  equation  (55)  further  reduces  to 


>(«,  t)  _  ^ 


./*"T 


bmm 


A  further  reduction  results  from  letting  s  =  ^-A-  » 

Zj 


At 

1554.l-i.d-i)  l  +  4^ 

Z 1 10(1  2  \  1  —  S  /  1  + 
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where  71  =  0,1,  2, - ,  and 

\!hc  (2  nL-  x)  <  t  <  'Jhc  ^  2  (n  +  1)  L  -  *  J  . 

For  the  two  positions  of  interest,  x  =  o,  and  x  =  L,  the  solutions  are  as 
follows: 

-  "+1  (58) 


for  x  =  o 


p(o,t)  _  1  d 
Z  \lV;  I 


2  n 


(1  +  S)n 

(wi(2("  +  1)  • 


for  x  =  L 


l  (i-S\n 

H\  V +  8) 

(virr  ((2“ 


(2  n  -  1) 


+  1) 


(59) 


correlation  of  lumped  and  distributed  parameter  analyses 

An  examination  of  the  distributed  parameter  frequency  response  discloses 
that  for  small  values  of  the  frequency  variable  «a  the  lumped  and  distributed 
parameter  responses  are  equivalent.  In  order  to  show  this,  equation  (27)  can 


be  written  as 


p  (/&>)  _  _ 

yR&W 


P_ 

W 


(/<d)  +  1 


(60) 


For  the  special  case  r  =  o  equation  (36)  can  be  written  as 


p(/<y)  A  +  jfi 


yi(/<y)  h  +  if* 


B 

Z  cos  Bx  +  / - sin  fix 

cos  B L  +  j  Z  sin  fiL 


For  values  of  <u  <<  — ,  cos/3x=l,  and  sinpx  =  px  because  Bx  1  /3L 
L 


o)  L 


«  1. 


Thus, 


7  ,  .  B  a  7  .  B*x 

p  (fa)  =  Z  +  >mPx  =  *  + 

«>((/<»)  I  ■  Z  coc  „L  1  +  j  Z  cocL 

'  B 


(61) 
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P  A  R 

Now,  remembering  that  /3  = —  for  r  =  o,  z  =  —  j  c  = — u2  =  ygRQ  and  LAn  =  V , 

u  w  u2 

equation  (61)  can  be  reduced  to 


p  (jco) 
u>i  (jco) 


P 


PV 

yR®W 


gx 

uC 


(jco)  + 


(62) 


1  P 

Since  /3*  « 1  and  — <  (Fig.  8),  equation  (62)  reduces  to  equation  (27). 

The  deviation  of  the  distributed  parameter  frequency  response  from  the 

lumped  parameter  begins  to  become  excessive  for  frequencies  just  below  the 

first  critical  frequency  (<u„=  -^—).  As  previously  noted  the  break  in  the  lumped 

0  2  L 

constant  frequency  response  occurs  at  =  —  =  -p®—  .  The  ratio  of  these  is 

yR®W 

m  ~W~  (63) 

coc  nu 

2L 

This  expression  reduces  to 

^b_  =2JLM  (64) 

_ *>C  _  _  _ 

where  M  is  the  Mach  number  of  the  gas  flowing  in  the  system.  The  constant  2yA 
is  about  unity  so  that  the  ratio  is  about  equal  to  the  Mach  number.  If  the  ratio  is 
small  the  deviation  occurs  far  down  on  the  frequency  response  gain  curve;  if  the 
ratio  is  large  the  deviation  occurs  early  on  the  frequency  response  plot. 

Examination  of  the  indicial  response  indicates  that  the  solution  (56)  can  be 
reduced  to  the  lumped  parameter  solution  (28)  by  the  following  treatment.  Sup¬ 
pose  the  ducting  length  L  tends  to  zero  while  the  volume  L  A0  remains  fixed. 

Then  tends  to  unity,  A0  tends  to  infinity,  V—  =  — —  tends  to  zero,  and  thus 

L  AoS 
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8  approaches  zero.  The  limit  on  time  may  then  be  written  as 

(2“- 

For  any  fixed  time  i  as  L  approaches  zero,  n  tends  to  — 


p  (*,  t)  A  -S\2Va7l 

z\»i\  v  +  ^  / 


(65) 


Then 


(66) 


Now,  remembering  that  Lim  (1  +  6*)*  =  e’° 

x  •*  o 


it  is  convenient  to  write 


8t 


p  (*,  t)  x  _  f  l  1  -  s\g  1 2  L 


-imt 


then 


ft#- 

Letting  6=1  and  6  =  -1  ,  respectively 


i  i 

~Tf  If 

(i  +  s)  a  -  s)  . 


(67) 


since  the  quantity  —  can  be  reduced  to  -  ,  which  is  the  time  constant 

..  8  yRBW 

t  ,  equation  (28).  Thus,  the  distributed  parameter  solution  can  be  reduced  to  the 
lumped  parameter  solution  if  the  length  of  the  ducting  is  made  small  and  the  area 
large. 
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OBSERVATIONS  ON  THE  RESPONSE  EQUATIONS 

There  are  a  number  of  interesting  observations  that  can  be  made  on  the  a- 
nalysis  so  far.  The  treatment  has  been  chiefly  analytical  without  regard  for 
physical  concepts.  This  section  is  devoted  to  observations  on  the  nature  of  the 
different  properties  of  gas  flow  systems  as  revealed  in  this  analysis. 

NATURE  OF  THE  FORCING  FUNCTION 

The  forcing  function,  or  disturbance,  considered  in  the  illustrative  example 
was  a  change  in  flow  into  the  system  at  the  inlet.  This  disturbance  is  the 
major  forcing  for  many  systems,  including  the  exhaust  system  of  an  engine 
test  plant.  One  might  as  easily  be  concerned  with  disturbances  in  inlet  or  out¬ 
let  pressure,  in  valve  area,  or  in  parameters  such  as  the  volume  of  a  bellows. 
In  any  case  the  method  of  treatment  would  be  similar  so  that  the  analysis  might 
be  used  on  such  varied  problems  as  pressure  response  of  instrument  lines  and 
transducers,  acoustical  design  of  mufflers,  analysis  of  pneumatic  control  ele¬ 
ments,  and  the  dynamic  analysis  of  wind  tunnel  performance.  The  specific 
discussion  is  limited  to  changes  in  flow  or  valve  area. 

NATURE  OF  THE  RESPONSE 

The  response  desired  of  gas  flow  systems  is  nearly  always  the  pressure 
change  due  to  some  forcing.  In  system  synthesis  for  control  purposes  the  steady- 
state  response  to  sinusoidal  forcing  and  the  time  response  to  step  function 
forcing  are  usually  sufficient  to  adequately  describe  control  characteristics. 
Depending  on  the  nature  of  the  problem,  varying  degrees  of  completeness  of  the 
response  may  be  needed.  For  some  problems  only  the  lumped  parameter  trans¬ 
fer  function  may  be  required.  For  others  the  complete  distributed  response  to  a 
step  or  other  function  may  be  required. 
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The  Lumped  Parameter  Transfer  Function 

In  the  description  of  systems,  perhaps  the  most  informative  knowledge  is  an 
expression  of  the  linear,  rion-dimensionalized  transfer  function.  One  sees  in  this 
expression  the  sensitivity  of  the  system  (gain),  the  degree  of  complexity,  or  or¬ 
der,  of  the  system,  the  magnitude  of  the  time  constants,  and  the  degree  of  sta¬ 
bility  to  be  expected.  Fairly  accurate  knowledge  of  these  is  in  many  instances 
sufficient  for  the  design  of  control  systems. 

Formulation  of  the  lumped  parameter  transfer  function  is  illustrated  by  the 
system  shown  in  Appendix  2.  The  equations  governing  lumped  parameter  are 
shown  at  their  respective  locations,  and  the  transfer  function  resulting  from  the 
combination  of  these  equations  is  also  shown.  In  the  system  of  Appendix  2,  a 
flow  of  gas  Wi  enters  the  volume  Vx,  passes  through  a  control  valve  system  into 
volume  V3l  and  is  pumped  out  by  an  exhauster  which  has  a  recirculation  line 
around  it.  The  control  valve  system  consists  of  a  throttle  valve  and  a  bleed 
valve.  The  inlet  flow  is  an  independent  quantity  which  deviates  from  a  steady 
state.  The  pressure  response  to  this  deviation  is  desired.  The  simultaneous 
solution  of  all  the  transformed  equations  which  describe  each  system  element 
yields  the  transfer  function. 

There  are  several  interesting  properties  of  this  transfer  function.  The 
system  contains  two  energy  storage  elements,  the  two  volumes;  and  hence 
a  second  order  denominated  in  the  transfer  function  is  assured.  The  system 
constants  are  all  real,  positive  numbers;  and  this  assures,  for  the  second 
order  system,  that  the  system  will  be  stable.  Furthermore,  when  the  trans¬ 
fer  function  is  written  in  the  form  shown,  where  a>n  is  the  natural  frequency 
and  C  isthe  damping  ratio,  itcan  be  shown  that  £  is  always  greater  than  unity, 
thus  assuring  that  the  system  is  over -damped.  This  also  assures  that  the 
denominator  has  two  real  roots  and  can  be  expressed  as  the  product  of  two  first 
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order  lags,  (r,  s  +  1)  (r4  s  +  1)  .  Now  it  can  also  be  shown  that  t3  <  r,  <  r4  and  that 
r,<  This  guarantees  that  the  log  magnitude  -  log  frequency  curve 

never  has  a  positive  slope  and  that  the  phase  shift  is  never  greater  than  -90°  for 
the  two  responses  to  valve  are  a  changes  and  the  response  to  flow  disturbance.  For 
the  surge  valve  area  change,  again  the  amplitude  frequency  response  always  has 
negative  slope,  but  the  phase  angle  reaches  -180°  at  infinite  frequency. 

Two  special  cases  which  occur  in  systems  of  this  kind  will  now  be  consid¬ 
ered.  In  Appendix  1  it  is  shown  that  for  pressure  ratios  across  a  restriction  less 
than  about  0.  528,  the  flow  is  no  longer  dependent  on  downstream  pressure;  that 
is,  in  Appendix  1,  kx  =  o  and  in  the  present  problem,  k,  =  o.  This  is  the  con¬ 
dition  of  a  "choked"  control  valve  in  the  present  system,  which  results  in  a  great 
simplification  of  the  transfer  function.  Under  this  condition  the  transfer  function 
consists  of  first  order  lags  as  shown  in  Appendix  2. 

A  second  special  case  which  is  often  encountered  occurs  when  the  constant 
ke  is  infinite.  One  problem  to  be  solved  in  gas  flow  systems  is  the  need  for  a 
safety  control  which  serves  to  keep  the  exhauster  away  from  the  surge  region. 
This  is  accomplished  by  the  recirculation  surge  valve  shown  in  Appendix  2. 
Suppose  this  valve  allows  sufficient  flow  so  that  the  exhauster  is  just  at  surge; 
that  is,  ke  =  oo .  Then  the  system  is  again  characterized  by  first  order  lags. 
(It  should  be  noted  that  the  problem  of  surge  control  can  be  attacked  using  the 
s  ame  governing  laws . ) 

The  Distributed  Parameter  Transfer  Function 

The  distributed  parameter  description  of  plant  response  adds  what  the  lumped 
parameter  description  fails  to  express- -the  spatial  effects.  Perhaps  the  most 
enlightening  picture  of  the  distributed  system  can  be  drawn  from  an  examination 
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of  the  physical  interpretation  of  the  solution  to  the  partial  differential  equations. 
Equations  (32)  and  (33)  and  more  generally,  equations  (45)  and  (46)  express  this 

general  solution.  Allowing  the  Laplace  operator  in  equation  (46)  to  take  on  its 

y 

general  complex  form,  the  term  (he  s 2  +  res)  3  can  be  expressed  as  the  com¬ 
plex  number  a  +  // 3  .  Then  equation  (46)  can  be  written  as 

.  P  (5)  =  E*  eax  eiPx+  F*  e~ax 

The  first  term  of  this  expression  represents  a  quantity  that  increases  in 
magnitude  and  increases  in  phase  as  x  increases.  This  term  must,  then,  rep¬ 
resent  the  contribution  to  the  pressure  change  which  originates  at  the  point  x  =  L  , 
since  it  is  logical  to  surmise  that  as  a  pressure  disturbance  travels  toward  *  =  o 
from  x  =  L  it  experiences  an  attenuation  and  a  time  lag.  By  the  same  argument 
it  can  be  concluded  that  the  latter  term  in  the  expression,  which  decreases  in 
magnitude  and  phase  as  *  increases,  must  originate  at  the  point,  *  =  o.  Thus, 
if  a  disturbance  such  as  the  flow  disturbance  previously  mentioned  originates  at 
x  =  o ,  then  all  terms  in  the  solution  containing  e~ax  e~iPx  represent  the  origi¬ 
nal  disturbance  and  successive  reflections  from  x  =  o,  whereas  the  terms  con¬ 
taining  eax  e  $x  must  represent  the  successive  reflections  from  the  point  x  =  L. 

A  number  of  special  cases  of  the  distributed  parameter  solutions  may  be 
treated,  each  of  which  adds  to  the  understanding  of  a  gas  flow  system.  Most  of 
these,  however,  are  treated  with  other  purposes  in  mind  than  control.  Examples 
are  theclosedand  open-end  organ  pipe,  which  in  our  interpretation  represent  prob¬ 
lems  of  infinite  and  zero  terminal  impedance  Z,  the  characterized  line (Z 
and  the  infinite  line  ( e~ax=  o)  as  treated  in  electrical  transmission-line  studies 
and  in  the  pneumatic  instrument  problem  in  which  a  long  instrument  line  is  ter¬ 
minated  in  an  instrument  volume.  Most  of  these  treatments  consider  only  the 
resistanceless  case  such  as  was  done  in  connection  with  the  simple  system  of 
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Fig.  5.  Space  will  not  be  given  here  to  these  special  cases  since  they  are  amply 
treated  in  textbooks  and  reports. 

Consideration  of  the  resistance  term  has  not  been  discussed.  The  effect 
of  including  the  parameter  can  be  understood  most  easily  by  numerical  solution 
of  the  system  represented  in  Fig.  5.  This  has  been  done  in  Appendix  3,  where 
the  calculated  response  curves  are  given.  It  is  seen  that  the  effect  of  the  re¬ 
sistance  is  to  smooth  out  the  peaks  and  valleys  in  the  frequency  response.  The 
numbers  used  in  the  calculations  were  chosen  for  ease  of  computation. 

NATURE  OF  THE  PARAMETERS 

The  properties  of  gas  flow  systems  are  the  pneumatic  capacitance,  inertance 
and  resistance  terms  (C,  c,  h,  andr),  and  the  perturbation  constants --the  partial 
derivatives  evaluated  at  points  denoted  by 'A..  These  properties  determine  the  sys¬ 
tem  parameters,  which  are  the  gains,  time  constants,  and  propagation  constants 
of  the  response  equations.  Considerable  information  can  be  gained  about  gas  flow 
systems  from  a  study  of  these  properties  and  parameters. 


Lumped  Parameter  High  Frequency  Approximation 

In  equation  2  7,  for  the  region  where  rs  (or  r/tu)  is  large  compared  to  unity, 

then 

P  (ia>)  =  1  . 

to  i  ( ja> )  jCat 


The  magnitude  is 
straight  line  of  slope, 
erating  point  ( Pss  and 


and  the  log  magnitude  vs  log  frequency  plot  is  simply  a 

V 


-1.  Now,  the  property  C  ( 


yR& 


)  does  not  change  with  op- 


Jf'ss  ),  and  thus  the  high  frequency  portion  of  the  frequency 
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response  remains  fixed.  The  system  of  Appendix  2  exhibits  this  same  property 

also,  where  it  can  be  shown  that  for  high  frequencies 

Pi  (jo)  =  1  . 

U  'i  ( jco )  /  C  lco 

In  fact,  it  can  be  shown  that  for  all  gas  flow  systems,  using  only  the  lumped  pa¬ 
rameter  description,  the  high  frequency  region  of  the  frequency  response  remains 
fixed  and  only  the  low  frequency  region  changes  with  operating  point. 

For  ease  of  plotting  the  frequency  response,  it  is  convenient  to  know  the 
value  of  the  point  on  the  high  frequency  asymptote.  The  frequency  at  which  the 
high  frequency  approximation  log  magnitude  vs  log  frequency  plot  crosses  unit 
magnitude  is  ml  =  -k-  =  .  Now,  since  0  is  absolute  temperature,  tu,  does 

not  shift  appreciably  with  temperature.  The.  volume  is,  then,  the  determining 
factor  in  the  location  of  <u,  . 

The  Perturbation  Constants 

The  use  of  the  theory  of  small  perturbations  in  system  analysis  is  becoming 

/ 

widespread.  Probably  the  largest  field  of  use  is  in  the  aerodynamics  of  air¬ 
frames  where  the  partial  derivatives,  or  perturbation  constants  when  evaluated 
at  a  point,  are  known  as  stability  derivatives.  The  evaluation  of  these  stability 
derivatives  of  airframe  dynamics  is  one  of  the  principal  projects  of  wind  tunnel 
testing.  The  perturbation  constants  of  the  valves  and  compressors,  for  accurate 
analysis,  must  also  be  taken  from  experiment.  However,  one  desires  to  be  able 
to  express  them  analytically.  For  a  number  of  systems  this  can  be  done,  as  is 
shown  for  a  control  valve  in  Appendix  1. 

The  Distributed  Parameters 

Pneumatic  resistance,  capacitance,  and  inertance  derive  their  names  from 

their  similarity  to  the  equivalent  electrical  terms.  The  form  of  the  resistance 

varies  greatly  with  the  system.  For  large  ducts  the  only  resistance  may  be  due 
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to  the  geometry  of  the  system,  such  as  tees,  elbows,  expansions  or  contractions. 
It  is  sometimes  sufficient  to  consider  the  total  resistance  as  being  distributed 
.  .  evenly  over  the  whole  length.  The  inertance  term  ( — - — )  becomes  much  larger 

4  e 

for  small  diameter  tubes  than  for  large  ducts,  while  for  the  capacitance  term 
A  £ 

( — o_£_)  just  the  opposite  is  true. 
u 

Through  the  entire  treatment  the  velocity  of  propagation  term  plays  a  major 
role.  This  is  true  also  for  the  lumped  pneumatic  capacitance  ( G  -  — —  =  —  - 

yR.©  u2 

'  and  thus  C  =  Lc)  .  The  equivalence  of  the  speed  of  sound  in  the  pneumatic  sys¬ 

tem  and  the  speed  of  light  in  the  electrical  system  should  also  be  noted. 

USE  OF  THE  ANALYSIS  IN  CONTROL  PROBLEMS 

A  number  of  factors  influence  the  manner  of  analysis  of  gas  flow  systems. 

* 

A  brief  cursory  analysis  of  a  system  which  already  exists  will  reveal  the  basic 
properties.  The  magnitude  of  the  largest  volume  and  rough  knowledge  of  tem¬ 
perature,  pressure,  and  flow  determine  the  largest  and  smallest  time  constants 
to  be  encountered.  The  shape  of  the  system  and  its  approximate  length  and  re¬ 
flection  times  quickly  determine  the  critical  frequencies;  and,  in  fact,  if  one 
knows  the  design  Mach  number  it  becomes  apparent  immediately  whether  other 
than  a  lumped  constant  analysis  need  be  considered. 

* 

In  the  design  of  new  systems  the  requirements  of  the  system  determine  what 
properties  should  be  designed  into  it  and  which  should  be  deleted  for  best  per¬ 
formance.  Thus,  if  it  is  required  that  the  set  point  pressure  in  a  system  change 

a 

rapidly,  the  time  constant  must  be  small  and  hence  the  volume  small,  whereas  for 
systems  in  which  it  is  required  that  pressure  changes  to  disturbances  be  small,  a 
large  volume  is  dictated.  In  all  cases  it  is  desired  that  the  volumes  be  short  in 
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length  and  large  in  cross  section,  thus  eliminating  distributed  parameter  effects 
as  much  as  possible. 

A  detailed  knowledge  of  the  system  transfer  function  with  fairly  accurate 
quantitative  knowledge  of  the  gains  and  time  constants  is  necessary  in  control  sys¬ 
tem  design.  If  /u>c  is  small,  the  lumped  parameter  expression  is  sufficient; 
if  u>b  /a ) c  is  large,  it  is  necessary  at  least  to  determine  a>c  .  Further  information 
about  the  distributed  parameter  response  is  interesting  academically  and  for 
understanding  gas  flow  systems  but  does  not  add  information  helpful  in  control 
system  design.  It  is  sufficient  to  know  that  at  a  frequency  equal  to  2 «c  the  phase 
shift  is  already  180°,  and  hence  a  stability  limit  is  reached.  Present  knowledge 
does  not  contribute  a  solution  to  this  stability  problem. 

A  convenient  and  informative  method  of  presenting  the  variation  of  plant  pa¬ 
rameters  is  to  "map"  them  on  a  pressure -weight  flow  plot.  Consider,  for  ex¬ 
ample,  the  time  constant  in  the  simple  system,  r  =  C/A;,.  This  maybe  reduced  to 

P  =  w 

which,  for  constant  temperature,  consists  of  a  family  of  straight  lines  on  the 
pressure-flow  map.  The  constants  of  the  system  of  Appendix  2  maybe  mapped 
similarly  although  the  relations  may  become  quite  complicated  at  times. 

With  a  knowledge  of  the  order  of  the  lumped  parameter  transfer  function, 
the  magnitude  of  the  gains  and  time  constants,  and  the  distributed  parameter  cri¬ 
tical  frequencies,  the  control  system  designer  can  determine  his  controller  com¬ 
ponents,  the  speed  of  response  required  of  valves  and  sensing  elements,  and  the 
limits  of  the  control.  Then,  a  fairly  accurate  prediction  of  errors  and  speeds 
of  response  to  different  forcings  can  be  made. 
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APPENDIX  1 

COMPRESSIBLE  FLOW  THROUGH  A  CONTROL  VALVE 

The  controlled  area  of  a  restriction,  here  referred  to  as  a  control  valve, 
is  nearly  always  the  final  control  element  in  a  gas  flow  system.  The  system  de¬ 
signer  must,  then,  have  a  fairly  accurate  knowledge  of  the  behavior  of  the  flow 
through  restrictions.  A  qualitative  knowledge  (e.g,,  the  flow  is  dependent  on  the 
pressure  ratio,  the  area  and  temperature)  is  not  sufficient;  it  is  also  necessary 
to  know  quantitatively  how  the  flow  is  dependent  on  the  variables. 

A  convenient  manner  of  expressing  this  dependence  is  by  a  measure  of  the 
deviation  of  the  actual,  measured  flow  from  a  general,  ideal  relationship.  That  is 

Va  =  Cf  Vt  ,  (68) 

where  Cf  is  a  "flow  coefficient"  and  subscripts  a  and  t  represent  actual  and  the¬ 
oretical  flows,  respectively. 

The  ideal  flow  relationship  is  derived  in  most  texts  on  engineering  thermo¬ 
dynamics  (Ref.  3)  from  consideration  of  the  continuity  equation,  the  state  and 
process  equations  of  an  ideal  gas,  and  the  general  energy  equation.  A  brief 
derivation  of  this  ideal  expression  is  given  in  this  appendix. 

The  continuity  equation  can  be  expressed  as 


Px  4  vx  =  P  1  A1  =  Pi  4  b  . w  ,  (69) 

where  p  is  density,  A  is  open  area  and  v  is  velocity.  The  equation  of  state, 
p  =  P/RQ,  and  the  expression  for  an  adiabatic  change  of  state,  (p,/p*  )  =  KPJP%  )y  > 


are  also  of  use.  The  general  energy  equation  can  be  written,  for  the  stations 
i  and  *  as  shown  on  Fig.  9,  on  a  "per  pound  of  fluid  flowing"  basis,  as 

2  2 

Pt  V ,  -  Px  Vx  +  JQ+  Work  =  ](EX  -  Et)  +  JLJ~  >  (70> 
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FIG.  9.  SCHEMATIC  OF  RESTRICTION 


where  Q  is  heat  added,  E  is  internal  energy  and  /  is  the  heat-work  conversion 
coefficient.  If  the  process  is  adiabatic,  then  there  is  no  net  flow  of  heat  or  work. 
The  internal  energy  is  given  by  Cv  0,  and  by  use  of  the  equation  of  state  and  the 
adiabatic  relationship  the  energy  equation  can  be  reduced  to 

y-1 

1 


vx-vt 


(71) 


Remembering  that  R  =  /  cv  (y-1)  equation  (71)  further  reduces  to 

y-1 

,  /  2 vrA  „  r  _  ip: 

Vx  =  1 


OSH-©’  ].-■ 


Now,  using  the  continuity  equation. 


W 


=  Px  vx  =  J 


2 

y  +  i" 

7 

T  2  gy 

11 

fiV  (i 

kV 

V  PS  ( 1 

o  \y 
_£_)  L 

Lfi(y-l). 

®i 

Lvv  v 

v 

r2®s  \f 

\) 

> 

further  simplification  yields 


W 


=  A 


2ey  (p*\ 

-dV©, 


R(y 


) 


(72) 


(73) 


(74) 


If  A  ,  is  large  compared  to  Ax,  the  denominator  of  expression  (74)  reduces  to  unity 
and  the  equation  simplifies  to 


w 

1  2  gy 

Pi 

"k 

1 _ 

it 

X 
X  + 

k- » 

_ l 

A  ~ 

Ax 

\jR(y-l) 

V@7 

V\~  / 

■V  J 

(75) 
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The  function  of  the  pressure  ratio 


io' 


,  when  plotted,  exhibits  a  maximum 

y 


which  can  be  found  by  differentiation  to  occur  at 


_  (  _2_V  -  1 

=  \Y  +  1/ 


(76) 


FIG.  10.  PLOT  OF  FUNCTION  OF  PRESSURE  RATIO 


This  value  of  the  pressure  ratio,  when  substituted  in  equation  (72)  yields 


and  since' 


’■  -8*  fe) 


Y-  1 


/ 2 ygR®  2 

=  / _ L_  +  v.1 

V  y  +  1 

=  ®x  (?+■ 1  j  and  vx  «vx  for  <4,»  ^.then 


(77) 


%  =  yjyeR9x  (78) 

which  is  the  sonic  velocity.  Once  sonic  velocity  is  achieved,  no  change  in  pres¬ 
sure  occurring  downstream  of  point  *  can  be  transmitted  upstream  and  hence 
flow  is  independent  of  pressure  ratios  smaller  than  that  indicated  in  equation  (76). 

The  flow  for  "choked"  flows  then  reduces  to 

W  k  P, 


V0, 


(79) 


It  is  reasonable  to  suppose  that  the  flow  through  a  restriction  will,  after  a 
fashion,  obey  the  relationships  (75)  and  (79).  In  fact,  it  seems  that  the  flow  re¬ 
lationship  of  an  actual  valve  should  differ  from  the  theoretical  expression  by  only 
a  constant,  the  flow  coefficient.  That  is,  from  equation  (68)  , 
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'■  "  c'k  ~W~  1  (£)  m 

where  k  and  f  (~p~j  ^ave  the  obvious  definitions. 

To  be  valid,  this  approach  to  the  problem  should  be  checked  by  experiment 

for  each  type  of  control  valve.  This  check  has  been  made  for  a  "butterfly"  type 

control  valve  with  very  good  agreement  (Ref.  2). 

Equation  (80)  may  be  written  as 

Wa  =  f{Ax,Qlt  Px>  P,_)  . 

Taking  the  partial  differential  of  Wa  we  obtain 


dW  ja  ^  dW  dW 

mT  ®‘  *  w; 


or,  by  making  certain  definitions  tv  =  ka  a  +  k$  6  -  kx 
where 


a  =  dAx 

e  =  d® 

Px  -  dPx 
Pi  =  dpi 

k„ 


dw  _  w 


dAx  4 


ke 


ew 

d® 


j_  JL 
2  0 


it,  ♦  Wl 

dpl  (81> 

Xpl 

(82) 

M- 

*r-(£)  ” 

K  = 


dW 
dP , 


(PA 

\t,i 

2 

|y  - 

T  1 

(px) 

\p  1  j 

y  +  1 

1  y 

\  y  ) 

im 

The  partial  derivatives  as  expressed  here  suggest  a  convenient  method  of  quanti¬ 
tative  evaluation  of  the  perturbation  constants  at  an  operating  point,  W  and  P,  . 
The  quantities  kx  and  can  be  conveniently  plotted  in  a  dimensionless  form  as 
kpl/W  vs.  px/Pi .  It  is  also  convenient  to  plot  the  ratio  of  weight  flow  to  choked 


weight  flow  y~  which  is 


(?) 


choked 


which  is  X  can  also  be  evaluated,  since 
A 
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The  functions  and  kx  are  plotted  in  nondime  nsionalized  form  against 

Px  /Pl  on  Figs.  11a,  lib  and  11c,  respectively.  It  is  of  interest  to  note  that  the 
flow  coefficient  is  absent  from  these  expressions. 

It  should  be  observed  that  the  pressure  Px  is  the  "throat"  pressure  and  that 
the  pressure  downstream  of  the  throat  maybe  greater  than  Px.  This  is  consid¬ 
ered  as  a  "recovery"  and  introduces  some  error  when  throat  and  downstream 
pressure  are  considered  equal. 
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APPENDIX  2 

DETERMINATION  OF  THE  TRANSFER  FUNCTION  OF  A 
TYPICAL  GAS  FLOW  SYSTEM 

The  system  to  be  considered  in  this  appendix  may  represent  a  number  of 
different  systems,  but  the  author  is  familiar  with  it  as  the  exhaust  side  of  an 
engine-test  plant.  Referring  to  the  symbolic  diagram.  Fig.  12,  a  quantity  of 
gas  Wi  enters  volume  1,  which  is  some  geometry  of  ducting.  A  quantity  of  gas 

/ 

also  enters  volume  1  through  the  control  bleed-in  valve,  the  pressure  Pt 
being  less  than  Pa.  The  total  flow  then  enters  volume  2  through  the  throttle 
control  valve  and  is  exhausted  through  a  compressor  to  a  pressure  Pa  .  A  part 
of  the  total  flow  IFS  is  recirculated  around  the  exhauster  and  back  into  volume 
/  2  for  purposes  of  compressor  surge  control.  The  problem  is  to  determine  the 

^  dynamic  relationship  of  the  pressure  in  volume  1  as  a  function  of  an  inlet  flow 

disturbance  and  disturbances  in  the  three  valve  areas.  The  methods  of  the 

a 

lumped  parameter  analysis  given  in  this  paper  will  be  used.  Only  the  governing 
equations  and  the  final  transfer  function  will  be  given. 
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Conditions :  0,  =  02 

Pa  and  0a  are  constant. 

Governing  Equations: 


P*  = 


Wb  « 


W 


—  (Wi  +  Ffc  -  wt) ;  Pt  =  — -i - (w;  +  u>b  -  wt ) 

'-‘is  ,  .  Lts 

PPa/4i,  /  P.  \ 

—  i  + *“*s 

KP,At 


t  ~ 


vi r 


f 


ir:)-’"’- 


KP\  +  —  Pi 


p,  =  ~ —  ivt  +  ws  -  we ) ;  Pi  =  — («»*  +  »*  -  u,e  ) 

(-*2  ®  Cj  s 


=  f(P*)*  we  ~  keP* 


= 


KPa4  /  P, 


Vg^T 


^  1  r.  )  ;  WS  *■  &saS  “  ^«Pa 


(84) 

(85) 

(86) 

(87) 

(88) 

(89) 


Transfer  Function:  Response  of  P,  to  changes  in  Wi  ,  Ab  ,  A,  ,  As  is: 


P.(s)  = 


Kj  (r,-  s  +  1) 

J_  S2+  -l-C  5  +  1 
Wn1  <U„ 


w:  + 


Kb  (ri  «  +  1) 


1  s’  +  JLt  s  +  1 


ab  _  **  (r»  S  +  1} - a, 

-i-52  +  JL(  s  +  1 


<*)„  <un 


- *i_ - „ 

-V  t  — —  (  s  +  1 


or , 


Pi  (s)  = 


K-i  (r,  s  +  1)  w  Kb  (’’Is  +  1  )  Kt  (r2  s  +  1)  Ks 

(r,  s  +1)  (r4  s +.1)  1  +  (rs  s  +  l)(r4  s  +  1)  b  (r3  s  +  l)(r<  s  +  1)  °(  +  (r,  s  +  l)(r4  s  +  1)  S 
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Parameters: 


Ki 


h  = 


K 


t  = 


K, 


ks  + 

ks  +  k 

e 

k2ks 

+  ktk5  + 

Me  + 

k2k. 

+  k3kg 

kb  (k,  - 

h  ks  + 

ke) 

k2k2 

+  ktks  + 

Me  + 

Ms 

+  Kkg 

kt  (  ks 

+  ke) 

kA 

+  ktks  + 

Me  + 

+  kjkg 

k$  k$ 

1 


v^c: 


V&J&3  +  fcjAg  +  AlAe  +  AaAj  +  A,Ae 


£  .  1.  V 


Ct 


(t5  +  &j  +  &_)  + 


+  K) 


■\J~icJc,  +  icjcs  +  kjc^+  k2ks  +  k2ke 


C, 


k,  +  k,  +  ke 


;  Ti 


k,  +  kg 


T.,rt  = 


”  '  ■ 

^5  +  +  &2 

tl/s 

(k,+k,+ke  kl  +  k2) 

4 

c,  c, 

\  c,  .  cx  ) 

+ 

C,C2 

Special  Case  of  a  Choked  Control  Valve: 
k,  =  0 


r,  =  r,  =  r,  ;  r4 


C. 


4,  +  &3 


if*  =  ■  -  1  ;  *6 - ^ - K,  =  - **- 

Aj  +  k2  +  k2  '  kt  ■$■  k2 


Ki 


«)£  + 


«&  - 


K, 


PiU)  = 

wi(s)  r4  s  +  1  r4  s  +  1  '  r4  s  +  1 


t _ a, 


,  h  s  =  0 
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Special  Case  of  Exhauster  at  Surge  Flow: 


1  2£  C, 

r,  =  r2  =  0;  -  =  0,  -2-  =  - — l—  =  r„ 


< » 


con  A.+A, 


K:  =  - 

A,  +  A3 


Ku  = 


A:,  +  Aj 


**  = 


A. 


A,  +  A, 


Pi  («) 
«J;  (s) 


Kj 

r,s  +  1 


*6 

+  - 

r,s  +  1 


ab 


K, 


r.s  +  1 


at 
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APPENDIX  3 


NUMERICAL  CALCULATION  OF  THE  RESPONSE  OF  A  SIMPLE  SYSTEM 


In  order  to  show  more  definitely  the  response  of  a  system  and  the  effect  of 
different  properties  of  the  system,  the  response  of  the  system  of  Fig.  13  is 
computed.  The  numbers  have  been  chosen  so  as  to  minimize  the  computations. 
The  ft-lb-sec  unit  system  is  used. 


Properties 

Duct  length,  L  =  50  nr  ft 
Duct  cross-sectional  area,  A0 


2000 

ng 


ft 


Steady  state  pressure,  Pss  =  1000  lbs/ft 2 
Steady  state  flow,  B®®  =  100  Ibs/sec 
Steady  state  temperature,  0®®  =  416°/? 

Velocity  of  wave  propagation,  u  =  1000  ft/sec 
Lumped  pneumatic  capacitance,  C  =  — =  0.1  ft2 

yfi  W 

Distributed  pneumatic  capacitance,  c  =  -^*2.  =  ~  x  10'*  ft 


sec1 


Distributed  pneumatic  inertance,  h  =  xy  =  -j— *  10"  IX- 
Distributed  pneumatic  resistance,  r  =  x  10" 


53 


AEDC-TR-55-11 


Parameters 

Terminal  pneumatic  impedance,  Z  =  ~  =  10  sec/ft2 

•  r  •  K-O  • 

Lumped  parameter  time  constant,  r  =  — =  1  sec 

Kq 

Distributed  parameter  critical  frequency,  coc  =-|jr-  =  10  rad/sec 

Lumped  Parameter  Frequency  Response 

The  lumped  parameter  transfer  function  can  be  written  numerically  as 

P  (s)  1  10  (90) 

Wi  (s)  k0  (rs  + 1)  s  +  1 

or  replacing  s  by  jco ,  the  frequency  response  is 

P  (M  10  (  91) 

w-  (jco)  jco  +  1 

The  frequency  response  is  shown  on  Figs.  14  and  15  for  the  log-log  plot  of 
magnitude  ratio  vs  frequency  and  the  semi-log  plot  of  phase  angle  vs  frequency, 
respectively. 


Distributed  Parameter  Frequency  Response 

The  distributed  parameter  frequency  response  transfer  function  can  be 
written  as 


p  (*,  jco) 
U>1  ( X ,  jco) 


h  +  ifx 
/»"+  iU 


(92) 


where  the  functions  i  are  given  by 

/,  =  Z  cos  (/3*)cosh  ( ax )  +  —  cos  (Bx)  sinh  (ax)  +  -  sin  (Bx)  cosh  (ax) 


u 


Z  sin  (Bx)  sinh  (ax)  +  — —  sin  (Bx)  cosh  (ax)  - 
coc 


—  cos  (Bx)  sinh  (ax) 

coc 


/,  =  cos  (BL)  cosh  (aL)  +  ^,^1  cos  (BL)  sinh  (aL)  -  8>n  (BL)  cosh  (aL) 

=  sin  (BL)  sinh  (aL)  +  sin  (BL)  cosh  (aL)  +  ^  cos  (j8L)  sinh  (aL)  . 
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The  numbers  a  and  8  are  given  by 


[  ct>A  + 


Va>2  h2  +  r2  3 


a  =-&%£- 
2(8 

Case  1;  r  =  0 

Under  this  special  case  the  functions  simplify  greatly;  /3  =  ~-  and  a  =  0,  and 
the  functions  f  become 


fi 


Z  COS  8*  “  10  cos 


OX 

1000 


A  = 


-  JL 


o>c 


sin  8*  “  -5-  sin 


/,  =  cos  81<  “  cos  0.05  no 

t  ^  6jC  .  ^)r  20  ■  n  nf 

f.  =  — g —  sin  8^  “  -  sin  0.05  »ra> , 

8  it 

Two  positions  of  *  are  of  interest,  *  =  0  and  *  ■=>  L .  Under  these  conditions  the 
response  curves  are  shown  in  Figs.  14  and  15. 

Case  2:  r  4  0 

For  this  more  general  case  the  response  is  more  difficult  to  calculate.  Qur 
choice  of  r  makes  the  calculations  simpler  however.  Since  r  was  chosen  such 
that  r  =  h,  the  expressions  reduce  to 

8  =  ~~  v -j—  w  (tt>  +  V  o)3  +  lF 

1  I  2tu 
2u  +  1 


Point  by  point  calculation  of  the  response  then  yields  the  curves  shown  in  Figs. 
14  and  15. 

Lumped  Parameter  Indicia!  Response 

The  solution  of  the  transfer  function  (90)  for  wi  (s)  =  a  unit  step  function, 
is 

p(t)  10(1  —  e“ *) '  (93) 

This  response  is  plotted  as  a  function  of  time  on  Fig.  16. 
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Distributed  Parameter  Indicial  Response 


The  indicial  response  for  the  distributed,  resistanceless  system,  is  given 
by  the  expression 


p(x,  t)  =  Z  I 


1  - 


i  a  - s)  [  i  *  (r^-f-)1  ]  [ttt]  "  }  <94) 


(2n  L  -  x)  ^-~L.  Ain  +  1>  L  -  x  n  .  0,  1,  2,  —  , 


where  S  is  defined  by  the  ratio  _  .  For  the  two  positions  of  interest,  *  =  0 

Z‘ 

and  x  =  L,  the  responses  for  the  unit  step  input  are 
at  *  =  0 


P(0.  t)  =  x  (1-  8)w  +  1  _„1  (20  -  n) 


n  +  1 


2n 


(M 


(1  +  S)n 

2  (n  +  1) 


(20  +  n)n 


(95) 


at  x  =  L 


2Jl  1  A  -  8  Y  !  Ao  -  nY 

•t t  \1  +  8)  \20  +  nj 

(( 


(96) 


(2 n  -  1)  (  (  (2 n  +  1) 

\ 

This  plot  also  appears  on  Fig.  16 


Conclusions 


From  these  calculated  responses  one  gets  an  insight  into  the  nature  of  gas 
flow  systems.  The  effect  of  resistance  on  the  frequency  response  should  be 
particularly  noted.  It  would  be  well  also  to  have  the  indicial  response  for  the 
case  of  a  finite  resistance  but  this  calculation  would  be  very  difficult. 

The  following  conclusions  and  observations  can  be  made  on  the  frequency 
response  of  this  system: 

1.  The  departure  of  the  distributed  parameter  response  from  the  lumped 
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parameter  response  becomes  appreciable  at  about  half  way  between  the 
lumped  parameter  break  frequency  and  the  frequency  at  which  the  first 
quarter  wave  length  occurs. 

2.  The  effect  of  resistance  on  the  distributed  parameter  response  is  to  de¬ 
crease  the  magnitude  and  round  off  the  peaks  of  the  periodic  maxima  and 
minima. 

p 

3.  The  impedance  (equal  to  the  magnitude  ratio  ,  which  is  seen  at  the 
inlet  of  the  system)  is  increased  by  an  amount  equal  to  the  total  line  re¬ 
sistance  (r  L) .  The  outlet  impedance  does  not  experience  this  effect.  This 
is  shown  by  the  positions  on  the  curves  on  Fig.  14  at  ^  =  0.1  radians  per 
second. 

The  following  conclusions  and  observations  can  be  made  on  the  indicial 
response: 

1.  At  the  outlet,  no  effect  of  a  disturbance  at  the  inlet  is  felt  until  a  time 
0.05 ff  sec,  which  is  the  time  required  for  the  transmission  of  a  pressure 
wave.  Thereafter  steps  are  0.1*  sec  apart,  the  time  required  for  a  pres¬ 
sure  wave  to  traverse  the  length  in  both  directions. 

2.  The  pattern  of  the  successive  steps  shows  remarkable  resemblance  to 
the  first  order  response. 
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FIG.  15.  FREQUENCY  RESPONSE  ~  PHASE  ANGLE 
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FIG.  16.  1NDICIAL  RESPONSE 
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